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ABSTRACT 

Paramagnetic alignment of suprathermally rotating grains is discussed in view of recent 
progress in understanding subtle processes taking place over grain surface. It is shown 
that in typical ISM conditions, grains with surfaces of amorphous H2O ice, defected sili- 
cate or polymeric carbonaceous material are likely to exhibit enhanced alignment, while 
those of aromatic carbonaceous material or graphite are not. The critical grain sizes 
and temperature for the Purcell's alignment are obtained and preferential alignment of 
large grains is established. 
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INTRODUCTION 



The paramagnetic mechanism of grain alignment (Davis & 
1— Greenstein 1951) is believed to be too weak to account 
for observations unless enhanced imaginary part of grain 
I magnetic susceptibility is assumed (Jones & Spitzer 1967, 
! Mat his 1986). Physically this weakness means that ther- 
, mal disorientation due to chaotic atomic bombardment over- 



in' 

whelms orientation effects of dielectric grains with normal 
*** — . paramagnetic susceptibility rotating in the Galactic mag- 
*TH netic field. To diminish the effect of thermal disorientation 
ijit was suggested in Purcell (1975, 1979) that grains rotate 
Q suprathermally, i.e. with energy orders of magnitude greater 
than any temperature in the system times k (the Boltzmann 
X/l constant). 

It may be shown that among various causes of 
suprathermal rotation those associated with corpuscular and 
radiative fluxes result in the Gold-type alignment (Lazar- 
ian 1994a, 1994b); its time-scale is orders of magnitude less 
than that of paramagnetic relaxation. Therefore, 'chemically 
driven' suprathermal rotation is required by the Purcell's 
mechanism. The estimates in Spitzer & McGlynn (1979) 
showed that if high rates of grain resurfacing are consid- 
ered, it still improves the alignment only marginally. Thus, 
straightforward attempts to explain observations through 
the Purcell's mechanism (Aanestad & Greenberg 1983) had 
rather limited success (Spitzer 1993). 

The major difficulties of the paramagnetic alignment 
can be summarised as follows. 

(i) The paramagnetic alignment is too weak if ordinary para- 
magnetism is involved and 'standard' values of magnetic 
field and gaseous density are used (Spitzer & McGlynn 1979, 
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Mathis 1986, Chlewicki & Greenberg 1990, Greenberg 1993, 
Spitzer 1993). 

(ii) Observations indicate that grains of different chemical 
composition exhibit different degree of alignment (Mathis 
1986, Whittet 1992) while this does not follow from the the- 
ory (see Jones & Spitzer 1967). 

(iii) Large grains produce more polarization than small ones, 
while the Davis-Greenstein mechanism tends to go the other 
way (Mathis 1986, Kim & Martin 1994, Mathis 1994). 

In our recent paper on the alignment of fractal grains 
(Lazarian 1994c), we were mostly concerned with the first 
difficulty mentioned above. This papers studies details of 
grain surface physics and chemistry to address the other 
two problems. Such a study would not be possible until quite 
recently as it is based on new results obtained in the field 
(Duley 1993, Duley & Williams 1993). 

The structure of the paper is as follows. In section 2 
we discuss the physics and chemistry of suprathermal ro- 
tation caused by H2 formation over grain surfaces, whereas 
the paramagnetic relaxation is discussed in section 3. The 
ways of reconciling theory and observations are discussed 
in section 4 and two tests for the theory are suggested in 
section 5. 



2 'CHEMICALLY DRIVEN' SUPRATHERMAL 
ROTATION 

Although the interest to grain chemistry in view of grain 
alignment can be tracked back to Purcell & Spitzer's paper 
of 1971, it was only in Purcell (1975, 1979) that the concept 
of suprathermal rotation of grains driven by H2 formation 
was introduced. These pioneering papers were based on the 
level of knowledge about grain chemistry achieved at that 
time. New discoveries in this rapidly growing area make us 
to update some of the ideas put forward by Purcell. This 
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update does not alter the core of the Purcell's concept, but 
shows that the 'chemically driven' suprathermal rotation is 
not universally applicable, e.g. to grains of different chemical 
composition.^ 

It is shown in Duley & Williams (1993) that H 2 
molecules are formed over aromatic carbonaceous surfaces 
in low-energy state as the PAH.2H complex shares all the 
energy of stabilisation of H2 (4.5 eV) among its degrees of 
freedom'. This prediction entails that such grains should not 
rotate suprathermally due to H2 formation. 

The situation is different for H atoms weakly adsorbed 
by water ice or polymeric hydrocarbon. The binding is weak, 
its energy is not expected to exceed O.f eV (Lee 1972, Leitch- 
Devli'n & Williams 1984, Zhang & Buch 1990, Buch & Zhang 
1991). H2 molecules formed on such sites stabilise into very 
high vibrational states (Duley & Williams 1993). A portion 
of vibrational energy (~ O.f — 0.2 eV) is likely to be trans- 
ferred to the center-of-mass translational motion (Tielens & 
Allamandola 1987) at approximately I0 r s -1 rate (Lucas & 
Ewing 1981). Thus an injection of high speed H2 molecules 
is expected from icy or polymeric hydrocarbon grain sur- 
faces. Similarly using results in Duley & Williams (1986) it 
is possible to show that H2 molecules formed over defected 
silicate surfaces are also expected to have high kinetic ener- 
gies. For the hydrogenate amorphous carbon (HAC) model 
of interstellar dust (Jones et al. 1990) which consists of a 
mixture of polymeric, aromatic and diamond-like regions of 
carbon separated by voids, the situation is more involved. 
H2 formation is expected to be efficient in the sense of recoils 
deposited with grains for polymeric and probably diamond- 
like regions while aromatic regions are expected to eject H2 
with low kinetic energy. 

The ejection of H2 molecules with high kinetic energy 
E results in torques applied to a grain. Everywhere below 
we assume E ~ 0.2 eV. This means that our approach is 
applicable to suprathermal rotation of grains with only icy, 
polymeric carbonaceous or silicate surfaces. 

It was shown in Purcell (1979) that internal dissipa- 
tion of energy within grains, mainly due to the Barnett 
relaxation, + suppresses rotation around any axis but the axis 
of the greatest inertia on the time-scale ~ 10 7 /?? s, where n 
is the ratio of grain rotational energy to the equipartition 
energy ~ kT. Thus grains rotate around their major axes of 
inertia that we denote here the 2-axis. Therefore the prob- 
lem becomes one dimensional as the two other components 
of the torque contribute only to insignificant nutations. The 
number of H2 molecules ejected per second from an indi- 
vidual site is ~ 71 I 2 vhuhv -1 , where 71 the portion of H 
atoms with velocity vh and concentration ran adsorbed by 
the grain, while v is the number of active sites over the 
grain surface and I is the characteristic grain size which co- 
incides with the diameter for a spherical grain. Then the 
mean square of a residual torque from an H2 molecule of 
mass rau2 is 

t We will consider rotation induced by H2 formation, although 
some other chemical reactions might also contribute to the grain 
rotation in dark clouds, where the abundance of atomic hydrogen 
is low. 

+ A quantitative treatment through solving Fokker-Planck equa- 
tions is given in Lazarian (1994a). 
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which spins up grains to angular velocities limited only by 
friction forces, i.e. 

fi = (M 2 } 1/2 ^ (2) 

where I z is the z component of the momentum of inertia, 
and td is the rotational dumping time (Spitzer & McGlynn 
1979): 

t d « QM gas = 0.6-^- (3) 
nmv\ 

where tgas is the time that takes a gain of density g s to 
collide with gaseous atoms of the net mass equal to that 
of the grain, ra and v\ denote the density and velocity of 
atoms, respectively, while m is the mass of an individual 
atom. To obtain both numerical values and functional de- 
pendencies further on we will use quantities normalised by 
their standard values. To tell the normalised values we will 
denote them by corresponding letters with hats, e.g. I is 
equivalent to 1/(2 • 10 -5 cm) with 2 • 10 -5 cm as a stan- 
dard value of a grain size. We consider 'standard' 71 = 0.2 
and VH2 = 5 • 10 5 cm s -1 . Speaking about diffuse clouds, 
we assume that all hydrogen there is in atomic form and 
therefore ra/rau = 1 and m/mu2 = 0.5. Note that some- 
times the choice of 'standard' values is arbitrary, e.g. for the 
time being we assume the density of active sites ctH2 to be 
10 11 cm -2 . However, we will see in section 3 that aH2 sensi- 
tively depends on subtle processes of poisoning and therefore 
'standard' values may be misleading. 

Using standard values of the parameters involved and 
assuming that the number of active sites is proportional to 
I 2 , it is possible to write the following expression for the 
angular velocity 

n « 1.4 • 10 8 71 a" 1 / 2 (fr 2 s (4) 

If grains are porous and fluffy as those found in the 
stratosphere (see fig. 7.3 in Williams 1993) the degree of 
suprathermal rotation can be limited provided that H2 for- 
mation takes place within narrow irregular pores of grains 
and nascent molecules thermalize due to numerous collisions 
with pore walls. These processes are discussed in more de- 
tail in Lazarian (1994c), and here we mention this in view 
of implications that the suprathermal rotation can have on 
grain structure. The tensile stress within the spinning grain 
is of the order of g s l 2 fl 2 may reach 10 8 dyn cm -2 . An in- 
teresting hypothesis by Wright (1994) is that suprathermal 
rotation limits the fractal dimension of grains may have the 
following consequences: 

a) if grains have high fractal dimension, their degree of 
suprathermality is low and their fractal dimension does not 
change, 

b) grains of aromatic carbonaceous material and graphite 
are not limited in their fractal dimension by suprathermal 
rotation, 

c) grains with surfaces of defected silicate, amorphous H2O 
ice or polymeric carbonaceous material with a limited num- 
ber of pores will have their fractal dimension decreased as a 
result of suprathermal rotation. 

In short, silicate, ice and carbonaceous grains may have 
either large d — > 3 or low d — > 2 fractal dimension. One 
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may also speculate that if ctH2 is proportional to grain sur- 
face, Q scales as l~ 2 [see Eq. (4)], the tensile stress scales 
as l~ 2 as well and therefore small grains are more vulner- 
able to centrifugal disruption, which may make them more 
spherical as compared to larger grains. Adopting a model 
of a grain as consisting of small particles jumbled together 
(Mathis 1990), one may say that these conglomerates are 
more subject to disruption if they are small. However, this 
problem needs further investigation and we will not proceed 
with its discussion here. 

Grain dynamics depends not only on the number of ac- 
tive sites but also on their distribution over grain surface. 
Both accretion of a new monomolyer and poisoning of ac- 
tive sites cause occasional 'crossovers' (Spitzer & McGlynn 
1979) which considerably effect grain alignment. If these 
crossovers occur on the time-scale /j,, one should multiply 
our estimates of Q given by Eq. (2) by t\ /2 / (t L + t d ) 1/2 . 



3 THE PURCELL'S MECHANISM 
3.1 Paramagnetic relaxation 

A paramagnetic grain rotating in the external magnetic field 
B is subjected to dissipative torques (Davis & Greenstein 
1951, see also Jones & Spitzer 1967, Purcell & Spitzer 1971, 
Roberge et al. 1993). In the grain coordinate system, the 
component of B perpendicular to the rotation axis appears 
as a rotating field. The magnetic moment produced by the 
imaginary part of magnetic susceptibility x" causes rota- 
tion of the vector of angular momentum towards the direc- 
tion of B. For low frequencies, x" °f paramagnetic materi- 
als is proportional to the component of the angular velocity 
perpendicular to B up to frequencies lo c ~ 10 — 10 8 s -1 . 
If 3> u>c, x" becomes independent of Q and approaches 
the static susceptibility. In short, large grain angular ve- 
locities increase the characteristic time of paramagnetic re- 
laxation as x" comes to its asymptotic values (compare 
Whittet 1992). The characteristic time of this relaxation is 
t r « lO^q^PfsB^Ki^e s (Spitzer 1978), where grain 
temperature T s and density q are normalised by 10 K and 
2 g cm -3 respectively, while magnetic field is normalised by 
2 • 10 -6 G, q\ is the enhancement factor of the imaginary 
part of magnetic susceptibility suggested in Jones & Spitzer 
(1967), and ii'(fi) is a function with the following asymp- 
totics 



1, 



(5) 



It was noted by Mathis (1994) that a low value of B is used 
above as standard and there might be a substantial increase 
for B 2 term (see e.g. Boulares & Cox 1990, fig 3b, where 
B 2 /8ir is shown in cgs units). 'However, many people would 
not be willing to go beyond about 3 micro G for the stan- 
dard ISM field in the Galactic plane', he concludes. Later in 
section 4 we will study the possibility that enhanced values 
of magnetic field are responsible for alignment. 

For simple estimates, one may consider that transition 
from one asymptotic to another is at Q = co c . Within this 
simple model, it is possible to introduce a critical grain size 
Icr m when the imaginary part of paramagnetic susceptibil- 



ity comes to saturation. For long-lived spin-up, i.e. /j, 3> td, 
Eqs. (4) and (5) provide 
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and for short-lived spin-up, i.e. ti <C td, it is (/i/Zd) 1 ^ 4 as 
large. 

For the original Davis-Greenstein process, t r should be 
compared with the time that takes a grain of mass m to 
collide with gas atoms of the same total mass tgas = 2 • 
10 12 Iq(vhAh)~ 1 s, where the gaseous density is normalised 
by 20 cm -2 and the 'standard' thermal velocity of atoms 
vh is assumed to be 1.3 • 10 5 cm s -1 . Note, that in many 
cases the velocity v is dominated by Alfvenic phenomena 
(see Aron & Max 1975, Elmegreen 1992, Myers 1987, 1994, 
Lazarian 1992, 1994a), the fact that makes tgas smaller. 

While a grain rotates suprathermally, chaotic bombard- 
ment of atoms cannot alter its rotation. It was argued in 
Purcell (1979) that a quantity t x = 1.3(/j, + td) should be 
compared with t r to find the Rayleigh reduction factor (a) 
introduced in Greenberg (1968), which relates the measures 
of polarization and grain alignment (see Whittet 1992). The 
expression of (a) is given in Aanestad & Greenberg (1983) 
and here we just mention that (a) depends on the following 
quantity (Spitzer & McGlynn 1979): 



(7) 



where F is the 'disorientation parameter' explicit analytical 
expression of which is obtained in the Appendix. One may 
establish a critical size l cr / after which F becomes greater 
than unity [see Eq. (All)] 




l cr f = 1.8 • 10" 
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,1/3 
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(8) 



where v is normalised by 100. If density of active sites is 
proportional to grain surface (we will see further on that 
this may not be the case): 



;i/ 4 



l cr f = 1.7 • 10" 



(9) 



For grains greater than this size, it is important to account 
for partial preservation of alignment during the crossovers. 

In the limit of /j, <C tgas and I < l cr /, 
8 e jj coincides with the parameter of alignment for the 
Davis-Greenstein mechanism and is equal to tgas jt r = 
0.02 (ifsVHnH)' 1 ^ 2 K(V)qi- 

Up to now we have discussed disorientation due to col- 
lisions with gaseous atoms. The disorientation due to radi- 
ation introduced in Purcell & Spitzer (1971) can be charac- 
terised by the time-scale t ra d ~ 5 • 10 13 lgT~ 4 s. 



3.2 Grain spin-up 

To obtain sufficient alignment for q\ = 1, it is required that 
th 3> tgas- The difficulty on this way is the finite life time 
of active sites over grain surface. 

A process that can limit the life time of active sites 
is resurfacing due to the accretion of a new monomolecular 
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layer (Spitzer & McGlynn 1979). A grain of mass m(l) grows 
as 

^± = ±S(l)vn a mS a (10) 

where f a is the adsorption probability, a coefficient with the 
range of values from to 1 for atoms of density n a , velocity 
v and mass m; S(l) ~ wl 2 is the geometric surface of the 
grain. Assuming that accreting atoms uniformly cover the 
entire physical surface of grains and substituting dm(l) = 
sS p h(l)dl, where S p h(l) is the physical surface area of grain, 
one obtains for the accreting time: 

f Ah\ fn H \ (m H \ f vh\ 1 
ta =lt gasq2 (—) (—) (—) (_)- (11) 

where Ah is the thickness of an adsorbed layer, and q 2 
is S p h(l)/ S(l) ratio. If vh and v are due to thermal mo- 
tions, then (■^-) = (7^") ^ an( i ta becomes proportional 
to (^f-) ^ • However, if non-thermal Alfvenic motions dom- 
inate, vh becomes equal to v and t a is proportional to (^f-) • 

We will see further on that oxygen is the most impor- 
tant agent poisoning the active sites. Other chemical el- 
ements either form compositions that desorb more easily 
than H2O or have a negligible abundance. As all heavy ele- 
ments, oxygen is usually depleted in the gaseous phase. For 
instance, its depletion towards £ Persei is -0.5, which gives 
wo /fiH ~ 10 -3 ' 7 (Cardelli et al. 1991). Assuming that a layer 
of Ah PS 3.7 • 10 -8 cm effectively covers active sites over 
grain surface (Aanestad & Greenberg 1983), one obtains 
Ah/l PS 1.85 • 10 -3 . Therefore for the accretion time-scale 
for oxygen is t a ss Stgas^a 1 ?2- The estimate for t a can be 
further increased if f a < 1 and the desorption is present. The 
value of f a ~ 0.1 was used in Spitzer (1978, p. 208) for illus- 
trative purposes while the values 0.02 and 0.2 were reported 
for sticking probabilities of C over silicate and graphite cov- 
ered with a monolayer of H2O surfaces (Williams 1993). It is 
shown in Jones & Williams (1984) that the effective sticking 
probability varies as mantles begin to be deposited and as 
a consequence the time of establishing monolayer coverage 
may be increased by a factor of about three. 

The products of grain surface chemistry are also sub- 
jected to desorption which increases the time of resurfacing. 
According to Jones & Williams (1984) the probability of re- 
taining the products OH and H2O on the grains surface is 
about 0.7 in dark quiescent clouds, whereas NH may be pre- 
dominantly released on formation (Wagenblast et al. 1993). 
Immediately after formation the reaction products are sub- 
jected to photodesorption. This process can be driven not 
only by UV quanta (Johnson 1982), but also by water ice 
3 /tin absorption (Williams et al. 1992). Note, that 3 /tin ra- 
diation penetrates deep into dark clouds. The energy release 
associated with H2 formation may be another cause of des- 
orbing weakly bound molecules such as CO in the vicinity 
of active sites (Duley & Williams 1993). This means that 
active sites can control to some extent resurfacing in their 
vicinity. In brief, we may expect to obtain t a of the order of 
10 tgas<l2 using rather conservative estimates for the above 
parameters. 

However, such considerations are applicable only if 
heavy atoms form an inert monolayer rather than poison- 
ing individual active sites. The thermal hopping of atoms 
has the characteristic time (Watson 1976) 



(12) 

where vq ps 10 12 s -1 is the characteristic frequency of vi- 
bration for an adsorbed particle and Eh is the energy of po- 
tential barrier which is approximately one third of the total 
binding energy (although see Jaycock & Parfitt 1986). The 
total binding energy for oxygen may be assumed ~ fc-800 K 
(Tielens & Allamandola 1988). In this case, the normalised 
E h is equal to E h /(k ■ 267 K). 

An adsorbed H atom scans the entire surface of a grain 
of I = 2-10 -5 cm on the time-scale not exceeding 10 -2 s. If it 
meets another H atom they are likely to form H2 molecule on 
a time-scale much shorter than the time of adsorbing another 
hydrogen atom. Therefore if an oxygen atom is adsorbed 
after this H2 formation took place, it is free to scan the 
grain surface unless a new H atom is adsorbed. 

Consider an oxygen atom arriving at the grain surface. 
The sites of chemical adsorption can be either all filled with 
H atoms or be empty due to a recent H2 formation. The 
number of empty sites is limited. The activation barrier of 
10 3 K and a width of 10 -8 cm is usually assumed for reac- 
tions between physically and chemically adsorbed hydrogen 
atoms (Tielens & Allamandola 1988). This gives the time- 
scale for the reaction of the order 1.3 -10 -6 s. The migration 
time-scale depends on the properties of the grain surface. 
Quantum tunneling corresponds to the time-scale of the or- 
der 10 -12 s, although estimates of the order 3-10 -9 s can be 
also found in the literature (Tielens & Allamandola 1988). 
As we will see further, the above difference is crucial. In- 
deed, assuming the time for migration 10 -12 s, one obtains 
that before reacting, a hydrogen atom has chances of visiting 
Per = 10 6 sites of chemical adsorption, while this estimate 
falls to just Per ~ 500 if the other value for the migration 
time-scale is assumed. In the former case a grain with high 
density of active sites will have not more than one empty 
active site of H2 formation even if a grain is completely cov- 
ered by active sites; in the latter case up to ~ vjv cr (if 
v > Vcr) active sites may be empty. These empty active 
sites are the prime targets of an O atom hopping over the 
surface. If O atom is trapped in an active site, it will stay 
there and will be eventually hydrogenated with the reaction 
products blocking the access to the site. 

If we assume that the number of sites of physical ad- 
sorption is N p h and among v active sites one is empty, the 
probability of an O atom not to fill it as a result of an indi- 
vidual hop is 1 — N~ h x ; this probability for mj, hops decreases 
to (1 - iV p - 1 ) m ". If m h > 1 and N ph > 1, this expression 
can be approximated: 

(l-^r^exp(-g). (13) 

The number of hops mj, before an O atom hydrogenated 
is proportional to the ratio of the time of adsorbing a hy- 
drogen atom t c = 1.5 • 10 3 (7iWij{iiji 2 ) -1 s to the time th 
given by Eq. (12). If a time-scale for hydrogen diffusion is 
~ 10 -12 s for a grain of 10 cm with high density of active 
sites hydrogen atom will scan the entire surface of the grain 
and will react with an oxygen atom. If the mobility of hy- 
drogen is lower, the probability to hydrogenate O atom over 
a grain will be suppressed and mj, will be increased by a fac- 
tor exp( — vjv cr ). Therefore for sufficiently high vjv CT ratio 
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the majority of O atoms will eventually fill the empty active 
sites and the fact that the number of vacancies per grain 
increases from one to vjv cr only contributes to the effect. 
The exponential increase in poisoning of active sites when 
their number exceed v cr is likely to decrease their number 
to Per ■ The characteristic time of poisoning is as follows 

tp „ 2vtc (—\ (^L) 1 (14) 

KvhJ \ n J 1 — exp( — rnhf-Nph) 

where vo is the velocity of O atoms with the number density 
no , the last factor is inversely proportional to the probabil- 
ity of O atom to fill an empty active site and the factor 2 
reflects that in half of the cases O atoms will find that all 
sites are filled. For small mh/N p h ratio, the exponent can be 
expanded into series and t p PS 2v jrah t a . In other words, for 
tp > ta, the number of active sites should not be less than 
rrih, which is for v jv cr < 1 is equal to 

m h = « 4 • IO^i/Ah)- 1 !^ 2 ex V {-E h /f s } (15) 

th 

Therefore a grain of 2 • 10 -5 cm should have not less than 2 • 
10 3 active sites for t p to be greater than t a . This means that 
the time-scale for hydrogen migration should be less than 
10 -9 s. If v = v C r = 500 as in example discussed above, ti 
becomes 0.25t a which corresponds to ~ 2.5tg a sq2- Although 
less than our estimate for the time of spin-up limited only 
by monomolecular layer accretion, this estimate provides a 
substantial increase of ti for irregular (g 2 3> 1) grains. The 
functional dependence of ti as a function of grain fractal 
dimension obtained in Lazarian (1994c) also stays unaltered 
as t p ~ t a - However, the processes discussed here depend not 
only on the relative abundance of oxygen as it was assumed 
in Lazarian (1994c) but also on the grain temperature T s 
and the energy of the potential barrier Eh- Both values can 
be measured and therefore long-lived spin-up can be tested. 

Consider, for instance, grains covered by H2O ice. Sim- 
ulations in Buch & Zhang (1991) indicate that active sites 
(i.e. sites with binding energy in excess of fc-900 K) cover 
around 10 -2 of the entire grain surface that a 'standard' 
grain will have v ~ 1.5 • 10 4 active sites. Data in table 2 in 
Tielens & Allamandola (1988) show that migration scales for 
hydrogen and oxygen are, respectively, 10 -12 s and 10 -2 s. 
For such high mobilities of hydrogen v cr > v and therefore 
v should be used in Eq. (14). A substitution of these values 
into Eq. (15) provides mj, PS 1.5 • 10 5 and thus ti becomes 
ten times less than t a - To obtain ti 3> td, one needs assume 
either }2 > 1 or T s < 10 K. For T s = 15 K, mj, becomes 
~ 6-10 7 and long-lived spin-up becomes improbable. On the 
basis of this data it is evident that the Purcell's alignment 
of grains with surfaces covered with amorphous H2O ice is 
limited to a rather narrow range of grain temperatures if our 
data for v and Eh are correct. 

The critical size of the grain corresponds to t c becom- 
ing equal to the time for an oxygen atom to scan the grain 
surface. The time t c varies as 7 -1 l -2 , the time for O atom 
to migrate over the grain surface as (j^ 2 ; hence the critical 
size is 

l crp w6-10- 6 (g 2 « ff 7 1 )-7f/exp {-0.25E h /f s } cm (16) 

If T s increases or Eh decreases, rapid poisoning becomes 
essential for grains of all sizes. At the first glance, l cr p 
seems to be insensitive to the atomic hydrogen concentra- 



tion: even a steady-state value of ran ~ 1cm -3 (see e.g. 
Duley & Williams 1984) the critical size becomes only twice 
as large as the one given by Eq. (16). However, this is only 
true if the atomic hydrogen concentration is not less than 
that of oxygen. Otherwise, there is too little of hydrogen to 
immobilise poisoning O atoms. 

For grains of size less than l cr p , the crossover time [see 
Eq. (14)] is in order of magnitude ti ~ vt c (nH /no). This 
time is ~ 10 8 for standard values of parameters adopted in 
the paper and is much less than the corresponding tgas ~ 
5 • 10 11 s. If the number of active sites is proportional to I 2 , 
than ti/tgas ~ 10 -3 will not change as the grain becomes 
smaller. Consequently, t x becomes equal to tgas and the im- 
provement of alignment as compared to Davis-Greenstein 
process is marginal. 



4 PROBLEMS AND SOLUTIONS 
4.1 Degree of alignment 

The most efficient polarising medium conceivable would con- 
tain perfectly aligned infinite cylinders with their axes per- 
pendicular to the line of sight. For such a model Mie cal- 
culations show that the ratio of polarization to extinction 
is of the order 0.3. At the same time, the observational up- 
per limit for the same ratio is < 0.064 (see Serkowski et al. 
1975). The difference between the two limits stems from the 
fact that realistic grains are irregular, moderately elongated 
particles and it is only a fraction of grains that is aligned. 
Nevertheless, it is difficult to account for the observations if 
8 e g does not exceed unity. 

Our estimates of 8 e g give ~ 0.3 for standard I and 
th = Wt gas (see section 2). Although 10 times greater than 
the initial value it is not sufficient to explain the polariza- 
tion. One of the ways to increase ti is by taking into ac- 
count photodesorption (see Johnson 1982). Indeed, as the 
characteristic time of photodesorption is around 10 9 s for 
unshielded regions of diffuse clouds, which is usually much 
less than the time of resurfacing, this may considerably in- 
crease th- However, the correlation of polarization with the 
strength of 3.1 /tin water ice feature (see Joyce & Simon 
1982) testify that the photodesorption cannot provide a uni- 
versal solution (see also Lee & Draine 1985). 

One may question 'standard' values for magnetic field 
that is responsible for grain alignment, while leaving the 
'standard' value of the mean magnetic field intact. Indeed, 
we observe the cumulative polarization and extinction along 
the line of sight, which presumably crosses quite different re- 
gions. If aligned dust grains are associated with clouds, we 
should use the data corresponding to the clouds. Zeeman 
measurements in Heiles (1989) provide B ss 1.2 • 10 -5 G 
for ra PS 7 cm -3 . This provides t r ~ 3 • 10 12 s, while 
tgas ~ 6-10 12 s, which gives for spherical grains tg a s/t r = 2. 
If one takes into account that grains are non-spherical (Pur- 
cell & Spitzer 1971) and subject to additional disorientation 
through emission, 8 e g ~ 0.8 can characterise alignment for 
th < td- This example of Heiles data for morphologically 
distinct shells may be considered to be extreme. However, 
numerous studies (see Myers 1985, 1987, Myers & Goodman 
1988) indicate that molecular clouds are magnetically sup- 
ported and exhibit motions which are supersonic 'in all but 
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the smallest cloud regions' (Myers 1985). These motions are 
usually attributed to Alfvenic turbulence (see Aron & Max 
1975, Falgarone & Puget 1986). However, the Alfvenic ve- 
locity is equal to the thermal velocity of atoms for 'standard' 
values of B = 2 • 10 -6 G and n = 20 cm -2 . Therefore to rec- 
oncile widely accepted explanations of the molecular cloud 
support and non-thermal line widths with observations, one 
should increase substantially B 2 j n ratio. Note that it is this 
ratio that enters the expression for tgas/tr and therefore an 
increase of Alfvenic velocity above the thermal one entails an 
increase of 6 e ff- H the Alfvenic velocity is five times greater 
than the thermal velocity 8 e g increases twenty five times as 
compared with 'standard' values. Even a less increase of 8 e g 
is adequate to explain the polarization if II is several times 
greater than tgas as it is discussed in section 3. 

Note, that the decrease of B 2 jn ratio as it follows from 
the data in Myers & Goodman (1988) may be one of the 
explanations for the decrease of polarization in dark clouds 
as indicated in Goodman (1994). 

Therefore the cumulative effects of inhibited poisoning 
as compared to that assumed Spitzer & McGlynn (1979) 
and enhanced values of magnetic field from recent observa- 
tions seems to favour the Purcell's mechanism to account 
at least for part of observational data. Indications that the 
degree of alignment correlates with the total magnetic field 
strength (Jones 1989), also favour such an interpretation as 
compared to the one where only grains with superparamag- 
netic inclusions are aligned and, if aligned, their alignment 
is perfect (Mathis 1986). 

The 'standard' values of the ratio B 2 jn appeared at 
the time when the Zeeman effect measurements were giving 
only upper limits of the magnetic field strength. According 
to Heiles (1989), this was an unfortunate result of the as- 
sumption that the splitting is easier to discover for narrow 
lines. In fact, strong magnetic fields usually correspond to 
broad lines. Present-day data both on Zeeman splitting and 
non-thermal linewidths tend to suggest that the 'standard' 
value of B 2 jn probably underestimate the actual value of 
this ratio. Having the long-lived spin-up with II within the 
range discussed in section 3, we feel quite comfortable with 
less than an order of magnitude increase of B 2 jn ratio as 
compared with its 'standard' value. 

Note that we are concerned with the ratio B 2 jn. For 
instance, using the 'standard' values of B = 3 • 10 -6 G but 
n = 5 cm -3 , one obtains t e g/t r K, 0.1, which is sufficient 
for explaining the polarization for the II /tgas ratio that we 
discussed. 

To summarise, the Purcell's mechanism seems to be ad- 
equately strong in the view of recent observational and the- 
oretical results. Its implications for grains of different sizes 
and chemical composition are discussed further on. 



4.2 Small and large grains 

Using the MNR grain model (Mathis et al. 1977), it is possi- 
ble to show (Mathis 1979) that a good fit of the polarization 
curve may be obtained if grains with radii above a certain 
threshold are aligned. The above dependence was accounted 
for in Mathis (1986) by assuming that only those grains 
that contain at least one superferromagnetic cluster are be- 
ing aligned. Within the adopted model, the probability of 



having such a cluster is greater for larger grains. Alterna- 
tively, for the Gold-type alignment (Lazarian 1994a), it is 
possible to show that only grains larger than 10 -5 cm are 
sufficiently inertial to be aligned under Alfvenic perturba- 
tions in the typical ISM conditions. § Here we are interested 
whether the Purcell's alignment can reproduce this type of 
behaviour for the polarization curve. 

It was noted in section 2 that if the density of active sites 
over the grain surface is low, the absolute majority of small 
grains would not have sites of H2 formation and therefore 
would not be aligned by the Purcell's mechanism. However, 
we do not believe that this is a probable explanation for the 
difference in the alignment of small and large grains. 

Our calculations in section 3 have shown that poisoning 
intensifies for I less than l cr p ~ 0.6 • 10 -5 cm. This critical 
size can further increase due to the inverse dependence of 
T s on I (Greenberg 1971): 

T s (l)xT s0 r 02 (17) 

This is a result of suppressed thermal emission from small 
grains in infrared. Although the power-law index is just 
0.2, this dependence may be essential in view of expo- 
nential dependence of l cr p on temperature [see Eq. (16)]. 
Thus, small grains are expected to be warm enough to have 
lcr p ~ 10 -5 cm. 

As soon as I < l cr p , rapid poisoning of active sites in- 
evitably causes II <C tgas - This forces down both the degree 
of suprathermality and the measure of alignment. One may 
erroneously conclude, that inspite of this drop, the degree of 
alignment has a formal limit of unity as I tends to zero. How- 
ever, t r ~ I 2 K(Q) and ii'(fi) ~ [see Eq. (5)] in the limit of 
high velocities, t r may be found to be inversely proportional 
to I if rotation is suprathermal and I < l cr p . At the same 
time, if poisoning of active sites for small grains is inhib- 
ited, then starting from l cr m t r becomes independent of I. 
Another important factor that suppresses the alignment of 
small grains is their disorientation due to thermal emission. 
In fact, one can introduce a critical size^ 

hr t ~ lO -6 ™^ 5 ^ 4 cm (18) 

starting from which the disorientation due to thermal emis- 
sion dominates. Due to the dependence of t ral j on and the 
change of grain temperature with I according to Eq. (17), 
the resultant dependence of t x PS t ral j is I 18 . Therefore for 

small g rams / <C lcr p and / <C lcr rrij $ eff ^ ^ 

At the same time, if rotation for small grains is not 
fast enough for their paramagnetic susceptibility to come to 
saturation, the degree of alignment increases with the de- 
crease of size as l~ ' 2 within the interval l cr m < I < lcr p- A 
transition from the long-lived spin-up to short-lived spin-up 
corresponds to a considerable drop in 8 e ff and therefore the 
alignment of grains with I < l cr p is expected to be marginal 
unless very high values of the ratio B 2 jn are present in the 

§ Grains are expected to be aligned with their long axes perpen- 
dicular to magnetic field lines, which corresponds to observations 
(Lazarian 1994a). 

^ This is not an exact result as grain emissivity deviates sub- 
stantially from the black body approximation for small grains. 
However, computations in Purcell &; Spitzer [1971] show that for 
small grains disorientation due to thermal emission dominates. 



Physics and chemistry of the Purcell's alignment 7 







J ^> Jcr m 


J <C Jcr m 


J ^> Jcr m 








J ^ Jcr t 


J ^ Jcr t 


J Jcr t 


J ^cr t 


/ J cr 


V 


,-2/3 


,4/3 


,-2/3 


,4/3 


J J cr 


f 










J J cr 


V 


r 2 


1° 


r 2 


1° 


J J Qy 


f 










J J Gf 


V 


,1/3 


jir/e 


,17/15 


,109/30 


J J cr 


f 










J J cr 


V 


r 1 


,3/2 


r i/ 5 


,23/10 


J J cv 


f 











Table 1. The dependence of 8 e jj on grain size for various combi- 
nations of the critical sizes if the number of active sites is propor- 
tional to the grain surface. l cr m = critical size below which the 
paramagnetic susceptibility is saturated [Eq. (6)], l cr $ = critical 
size above which a partial preservation of alignment in the course 
of crossovers is important [Eq. (9)], l cr p = critical size below 
which the oxygen poisoning of active sites is enhanced [Eq. (16)] 
l cr t = critical size below which disorientation due to thermal 
emission dominates [see Eq. (18)]. 
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Table 2. The dependence of S e jj on grain size when the number 
of active sites is stabilised at the level v CT . 

grain environment. Both the disorientation through radia- 
tion and saturation of paramagnetic dissipation suppress the 
alignment for the smallest grains. Our computations show 
that we should expect the said drop for grains less than 
10" 5 cm. 

For large grains, partial preservation of orientation dur- 
ing crossovers becomes important and this improves their 
alignment: for I > l cr m and I < l cr m , 8 e g becomes propor- 
tional to I -2 / 3 and I 4 / 3 , respectively. Other situations that 
may emerge for various combinations of four critical sizes 
and v < Per are shown in Table 1, whereas the dependen- 
cies for v > Per are given in Table 2. This information may 
be used for future quantitative comparison of theory and 
observations. 

Although possible in some situations, not all the combi- 
nations presented in Tables 1 and 2 are equally probable for 
diffuse clouds. We believe that it is more likely that there 



8 e g starts from ~ for the smallest grains and so as I ' 
until I = l C r m, after which it changes to I -1 / 5 A consider- 
able increase of 8 e g follows at I = l cr p , which corresponds 
to the transition to a long-lived spin-up. A decrease of 8 e g , 
which is proportional to I -2 , at I = l cr / becomes less rapid 
(~ I -2 / 3 ) due to a partial preservation of alignment during 
crossovers. When the number of active sites exceeds v cr and 
the decrease of 8 e g stops. 

It may be considered strange, but in view of new obser- 
vational results concerning magnetic fields in diffuse clouds 
Heiles (1989), Myers & Goodman (1988), Heiles et al. (1993) 
it may become more important to explain marginal align- 
ment of small grains rather than to account for good align- 
ment of large grains. Indeed, both the Purcell's alignment 
as it is discussed here and superparamagnetic alignment as 
it was suggested in Mathis (1986) fail to explain marginal 
alignment of smaller grains if the ratio B 2 /n is large enough. 

To summarise, the Purcell's mechanism is able to repro- 
duce qualitatively the observed features of the polarization 
curve. A more quantitative study will be done elsewhere. 

4.3 Grains of different chemical composition 

There exist a considerable convergence in views on the di- 
electric nature of polarising particles (see e.g. Whittet 1993), 
although the values of the imaginary part of the grain in- 
dex of refraction is a subject of controversy. Our results in- 
dicate, that 'H2 driven' suprathermal rotation and there- 
fore the Purcell's alignment are efficient for grains with 
surfaces of silicate, water ice or polymeric hydrocarbon. 
As the IR polarimetry shows substantial polarization in 
both silicate 9.7 fiia feature and 3.7 /tin water ice band 
(e.g. Whittet 1992) while graphite and amorphous carbon 
grains are likely not to be aligned (Mathis 1986), this may 
serve as an argument in favour of the Purcell's mechanism. 
Grains with surfaces of aromatic carbonaceous material and 
graphite should rotate thermally and therefore be subjected 
to the Davis-Greenstein alignment. But with the known ra- 
tios B 2 /n this alignment is marginal. If it is suprathermal 
rotation that controls the structure of grains (Wright 1994), 
the above grains are expected to be fluffy and therefore will 
not be aligned by Alfven-Gold processes discussed in Lazar- 
ian (1994a). 



5 TEST CASES 

5.1 Cold and warm grains 

We have seen that due to temperature increase, poisoning 
atoms become sufficiently mobile to attack active sites di- 
rectly and therefore the life time of active sites decreases. 
Our results show that for icy grains an increase of tem- 
perature up to 20 K leads to such a considerable increase 
in poisoning that long-lived spin-up becomes impossible. If 
the temperature increases further, H2 formation becomes 
suppressed. Indeed, for a reaction between physically and 
chemically adsorbed H atoms to occur, an activation bar- 
rier should be overcome. A critical temperature T cr can 
be defined for which the probability of reaction with a 
chemisorbed atoms is equal to the rate of evaporation: 
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k m\T r vo ) 

where r r is the reaction time-scale ~ 10 -6 s. Estimates pro- 
vide T cr ~ 50 K for H2 formation over a silicate surface 
(Tielens & Allamandola 1987). Above the critical tempera- 
ture, the reaction rate drops exponentially and the reduction 
factor at T s ~ 75 K is expected to be ~ 5 • 10 -3 . Such a de- 
crease in efficiency of H2 formation entails the same decrease 
in attainable angular velocity Q making it comparable with 
the velocity of thermal rotation lot- Therefore the Purcell's 
mechanism becomes inapplicable. Our discussion in section 3 
shows that even smaller increase of temperature makes long- 
lived spin-up impossible for grains with mantles of amor- 
phous ice. This means that, for instance, the aligned warm 
(T s ~ 75 K) dust in the core regions of OMC-1 (Chrysos- 
tomou et al. 1994) can not rotate suprathermally due to 
H2 formation. In general, it is unlikely that the suprather- 
mal paramagnetic alignment of warm (T s > 50 K) dust may 
be efficient. It is either the Gold-type alignment (see Lazar- 
ian 1994a) or Davis-Greenstein mechanism that governs the 
alignment of warm dust. To distinguish between the two 
cases one needs to study the polarization curve. The Davis- 
Greenstein mechanism should result in efficient alignment 
of small grains, II while Gold-type processes preferably align 
larger grains as a result of Alfvenic phenomena. 

5.2 H-atom population 

It is easy to see from Eqs. (1) — (3) that Q is proportional to 
ran/ra ratio and therefore becomes of the order of thermal 
angular velocity as soon as ran/ra ratio becomes as small 
as 10 -3 . Such circumstances are expected in dark clouds, 
where the UV is totally excluded and the residual H-atom 
population is maintained by cosmic rays; the steady state 
value of ran is ~ 1 cm -3 (Duley & Williams 1984). How- 
ever, there are indications (Norman & Silk 1980, Goldsmith 
et al. 1986) that this steady state may not be reached in 
many of dark clouds due to gas cycling between low and 
high density phases. For instance, ran in the range of 50- 
160 cm -3 is reported in Snell (1981) for several typical dark 
clouds which provides an estimate for ran/ra in the range 
of 2 • 10" 2 - 5 • 10" 4 , which can be interpreted as the evi- 
dence in favour of suprathermal rotation of grains in some of 
dark clouds. In general, suprathermal rotation and therefore 
the Purcell's alignment is expected to be more efficient in 
outer part of dark clouds, where ran/ra is greater. This qual- 
itatively corresponds to the results in Goodman (1994). In 
fact, we believe that both the decrease of B 2 jn ratio as a re- 
sult of ambipolar diffusion as well as ran/ra ratio contribute 
to the decrease of alignment in dark clouds. 



6 CONCLUSIONS 

A range of applicability of the Purcell's mechanism has been 
defined in view of new results in the field of grain physics 

II This is true unless the grains are Mathis' grains (1986), which 
population of larger grains is characterised by greater magnetic 
susceptibility. We will discuss this interesting proposal in our next 
paper. 



and chemistry. It has been shown that grains with amor- 
phous H2O ice mantles, defected silicate and polymeric car- 
bonaceous surfaces are likely to rotate suprathermally and 
can be aligned by the said mechanism if the relative abun- 
dance of atomic hydrogen ran/ra is not less than 10 -3 and 
grains are cool (T s < 50 K). At the same time, grains cov- 
ered with aromatic carbonaceous material or graphite are 
not expected to rotate suprathermally due to H2 formation 
and therefore this mechanism is not applicable to them. It 
has been shown that the Purcell's mechanism provides a 
preferential alignment of large grains, while only a marginal 
alignment of grains with I < 10 -5 cm is expected. 
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APPENDIX Al: 
PARAMETER 



DISORIENTATION 



A partial preservation of alignment during crossovers can 
be described using parameter F introduced in Spitzer & 
McGlynn (1979) as cos(xt) = exp( — F), where \T is the 
angle of deviation of the angular. Calculations in Spitzer & 
McGlynn (1979) give 



((AJ Z 



1 + 



3<(AJ ± ) 2 



2<(AJ Z 



(Al) 



F = - 

4 [\J L \(AJ Z 

in the approximation of the constant angular momentum 
component that is perpendicular to the axis of the greatest 
inertia. Note, that A J z and A Jj_ are elementary changes of 
the angular momentum in an individual torque event, and 
Jj_ is the value of the component of angular momentum 
perpendicular to the axis of major inertia. 

The squared dispersion of Jj_ at t = is given in Spitzer 
& McGlynn (1979, eq. 42): 

<Jl(0)> « {^f'\(^){A a Nfl 3 {AJ- 2 l 3 )ll 



T 2/3 



{(AJ ± 



(A2) 



where T(x) is the Gamma function, and A a is the coefficient 
of the Barnett relaxation (Purcell 1979, eq. 46) and N PS 
~fl 2 n,HVH is the frequency of the torque events. 

To calculate F from Eq. (Al), one needs to know the 
average of J^ J (0). Note, that ( Ji(0)) is a sum of (J 2 (0)} and 
(Jy(0)) and due to the symmetry inherent to the problem 
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For a Gaussian distribution, one gets 
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which can be calculated in polar coordinates to give 
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Similarly the Gaussian averages can be found for (A J z 2 ^ 3 ). 
The dispersion of (A J z ) which we will denote (72 is inversely 
proportional to the square root of the number of active sites 
v over the grain surface and depends on grain geometry and 
distribution of active sites in respect to the 2-axis. For a disc 
shaped grain (72 = 0.25lrriH2VH2V~ 1 ^ 2 ■ Therefore 
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(see Gradstein 
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[3.462(1)] and 

(A7) 



where D p (x) is the parabolic cylinder function and T(x) is 
the Gamma function. 

Thus substituting Eq. (A7) into Eq. (A2) 

<(AJ,) 2 } 



F 
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where the numerical factor is equal ti 




(A8) 




(A9) 



Assuming ((AJ Z ) 2 } = ((AJj_) 2 }, one obtains 

F= 5 -Z « A ff> 1/2 (A10) 

Substituting ((AJ Z ) 2 } = 15- "^2 "If 2 as wel l as expressions 
for A a , N and ([AJ Z ] 2 } into Eq. (A10), one gets 



where v is normalised by 100. If the number of active sites v 
over the grain surface is assumed to be proportional to the 
surface area of the grain S g {l), F decreases as I -4 / 3 with the 
increase of grain size. If the number of active sites is sta- 
bilised due to active poisoning at the level v cr , F decreases 
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